Extraction of level density and 7 strength 
q: -, function from primary 7 spectra 

A. Schiller* L. Bergholt, M. Guttormsen, E. Melby, J. Rekstad, 

and S. Siem 

Department of Physics, University of Oslo, 

P.O.Box 1048, Blindern, N-0316 Oslo, Norway 



Abstract 
> 

^\ • We present a new iterative procedure to extract the level density and 

f*"*) ' the 7 strength function from primary 7 spectra for energies close up to 

^^ . the neutron binding energy. The procedure is tested on simulated spectra 

and on data from the 173 Yb( 3 He,a) 172 Yb reaction. 
, ~~ l ' PACS number(s): 29.85. +c, 21.10.Ma, 25.55.Hp, 27.70. +q 
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X 1 Introduction 

The 7 transitions of excited nuclei give rich information on nuclear properties. 
In particular, the energy distribution of the first emitted 7 rays from a given 
excitation energy reveals information on the level density at the excitation en- 
ergy to which the nucleus decays, and the 7 strength function at the difference 
of those two energies. If the initial and final excitation energy belong to the 
^^ . continuum energy region, typically above 4 MeV of excitation energy for nuclei 

$H ' in the rare earth region, also thermodynamical properties may be investigated 

ii- 

Recently, the nuclear level density has become the object of new interest. 

There is strong theoretical progress in making calculations applicable to higher 
energies and heavier nuclei. In particular, the shell model Monte Carlo technique 
H, Q moves frontiers at present, and it is now mandatory to compare these 
calculations with experiments. Furthermore, the level density is essential for 
the understanding of the nucleon synthesis in stars. The level densities are 
input in large computer codes where thousands of cross sections are estimated 

Our present knowledge of the gross properties of the 7 strength function is 
also poor. The Weisskopf estimate which is based on single particle transitions, 
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see e.g. || , gives a first estimation for the strengths. However, for some measured 
7 transitions the transition rate may deviate many orders of magnitude from 
these estimates. A recent compilation on average 7 transition strengths for Ml, 
El and E2 transitions is given in Ref. H . The uncertainties concern the absolute 
strength as well as how the strength depends on the 7 transition energy. For El 
transitions, it is usually assumed that the energy dependency follows the GDR 
(7,7') cross section, however, this is not at all clear for low energy 7 rays. 

In this work we describe a method to extract simultaneously the level density 
and 7 strength function in the continuum energy region for low spin (0-6 h). 
The basic ideas and the assumptions behind the method were first presented in 
Ref. ||]. An implementation using an iterative projection technique, was first de- 
scribed in Ref. M. However, due to the existence of infinitely many solutions and 
the unfortunate renormalization of the primary 7 spectrum in every iteration 
step, this first implementation suffered from various severe problems, including 
divergence of the extracted quantities (To| . Several solutions of the convergence 
problem have been proposed and presented at different conferences, using ap- 
proximate normalizations, but none of them yielding exact reproductions of test 
spectra. However, data using one of those approximate methods were published 
in Ref. pi. Today, we consider the previous methods as premature, and we will 
present in the following a completely new, exact and convergent technique to 
extract level density and 7 strength function from primary 7 spectra. 

2 Extracting level density and 7 strength func- 
tion 

2.1 Ansatz 

We take the experimental primary 7 matrix T(Ei, E 7 ) as the starting point for 
this discussion. We assume that this matrix is normalized for every excitation 
energy bin Ei . This is done by letting the sum of T over all 7 energies E 1 from 
some minimum 7 energy -E™ ln to the maximum 7 energy Ei at this excitation 
energy bin be unity, i.e. 

Ei 

£ T(E l ,E 1 ) = l. (1) 

The 7 decay probability from the excitation energy Ei to Ef by a 7 ray 
with energy E-y = Ei — Ef in the continuum energy region is proportional to 



the level density g{Ef) and a 7 energy dependent factor F{E 1 ) (llj |12| , This 
ansatz is illustrated in Fig. |l|. The experimental normalized primary 7 matrix 



r can therefore theoretically be approximated by 



FjRy) g(E t - E 7 ) 






which also fulfills Eq. (|l|). 

As it is shown in Appendix [M one can construct all solutions of Eq. (0) by 
applying the transformation given by Eq. (H) to one arbitrary solution, where 
the generators of the transformation A, B and a can be chosen freely. 



g(Ei-E^) = giEi-E^AexpiaiEi-E^)), (3) 

F(E 7 ) = F(E y )Bexp(aE y ), 

2.2 Method 

2.2.1 th order estimate 

Since all possible solutions of Eq. (||) can be obtained by the transformation 
given by Eq. (H) of one arbitrary solution, we choose conveniently g(°> = 1. 
With this choice, the th order estimate of F is given by 



F^{E 1 



T(Ei,E^ = ^ Wi K 2L,„, - ( 4 ) 



Summing over the excitation energy interval Ef 1 ™ . . . £™ ax while obeying Ei > 
E 1 yields 

J2 r^E^^F^iE,) £ ' 



B i =max(£ min ,_B T ) £ l =max(_E min ,B 7 ) ^B 7 =-E™ ^ t ' 

(5) 
where the sum on the right hand side can be set to unity, giving 

£ max 

F(°)(25 7 )= £ r(^,E 7 ). (6) 

B i =max(_B i min ,B 7 ) 

2.2.2 Higher order estimates 

In order to calculate higher order estimates of the g and F functions, we devel- 
oped a least \ 2 method. The basic idea of this method is to minimize 

2 _ i y^ v^ / r th (^,g 7 ) -r(£ , t, J E 7 ) \ 

X N ficc ^ ^ \ ATiEuE-) J ' l j 



where N{ ree is the number of degrees of freedom, and AT(Ei, 2£ 7 ) is the uncer- 
tainty in the primary 7 matrix. Since we assume every point of the g and F 
functions as independent variables, we calculate iVf rGC as 

N ircc = ch(r) - ch(g) - ch(F), (8) 

where ch indicates the number of data points in the respective spectra. 
We minimize the reduced \ 2 by letting all derivatives 
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:X 2 -0 and — - — X 2 = (9) 



dF(E^)^ * dg(Ei-E^y 

for every argument £L and Ei — E^ respectively. A rather tedious but straight 
forward calculation yields equivalence between Eqs. (||) and 

T JE ] =maAEr ^ Ei) Q 2 (E l -E 1 )^{E^E 1 ) 

= z^ rax(£r ., g/+£rn) 1 /m {i (11) 
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where 
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s(Ei) -- 
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F(E y ) g(Ei - -E 7 ). 
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Within one iteration, we first calculate the functions a(Ei), b(Ei) and s(Ei), 
using the previous order estimates for g and F. Using these three functions, we 
can calculate the matrices (p(Ei,E^) and ip(Ei, E y ). Further on, we calculate 
the actual order estimates of g and F by means of Eqs. (Ilfl) and (O). Figure 
shows where the sums in Eqs. (|l0|) and (O) are performed. 



2.2.3 Convergence properties 

The method usually converges very well. However, in some cases the \ 2 mini- 
mum is very shallow, and the chance exists, that the iteration procedure might 
fail. In order to enhance convergence of the method, we have restricted the max- 
imum change of every data point in g and F within one iteration to a certain 
percentage P. This means that the data point obtained in the actual iteration 
(new) is checked if it lies within the interval 

< new < (1 + P/100) • old, (17) 



(l + P/100) 



determined by the data point from the previous iteration (old) . In the case that 
the new data point lies outside this interval, it will be set to the value of the 
closest boundary. 

Applying this method to some of our data, we have observed, that the smaller 
P is chosen, the smaller x 2 gets in the end, when the procedure has reached 
its limit. The reason for this is, that more and more data points in g and F 
will converge, while fewer and fewer points (typically at high energies £" 7 and 
Ef where few counts are available) are oscillating between the two boundaries 
given by Eq. (|l7|). Occasionally, we can choose P so small that all data points 
will converge and no oscillating behavior can be seen. However, in some cases 
oscillating data points can not be avoided by any choice of P which might 
indicate that the x 2 minimum is too shallow, or does not even exist, for some 
data points in g and F. 

A small P would lead to an accurate result but make a large number of 
iterations necessary, and a large P would shorten the execution time of the pro- 
cedure but affect the accurateness of the solution. We combine the advantages 
and avoid the disadvantages of the two concepts by letting P become smaller 
as a function of the number of iterations. In our actual computer code |L3], we 
have implemented a stepwise decrease of P as shown in Table (jL The choices 
of P as a function of the number of iterations is quite arbitrary, but we have 
achieved very good convergence for those spectra, where convergence properties 
without restrictions are rather fair. 

In conclusion we have to stress, that the convergence properties of the 
method in many cases do not require any restrictions of the maximum variation 
of data points within one iteration. In those cases however, where restrictions 
are mandatory to achieve or enhance convergence, they will only affect a small 
percentage of the data points at high energies, where data in the primary 7 
matrix are sparse and mainly erratically scattered. In those cases, where the 
restrictions of Table [l| would prove not to be satisfactory for convergence, the 
number of iterations or the value of P can be changed, since the validity of the 
method does not rely on these values. 



2.2.4 Error calculation 

A huge effort has been made in order to estimate errors of the data points in g 
and F. Since the experimental primary 7 matrix has been obtained from raw 
data by applying an unfolding procedure p4| and a subtraction technique fll5[ , 
error propagation through these methods is very tedious and has never been 
performed. In order to perform an error estimation of g and F, we first have to 
estimate the error of the primary 7 matrix data. A rough estimation yields 



Ar = 2 y/(M! + M 2 ) r, (18) 

where M\ denotes the number of first and higher generation 7 rays, and M 2 the 
number of second and higher generation 7 rays at one excitation energy bin Ei. 
We estimate those quantities roughly by 

Mi = max(l, M(Ei)) and M 2 = max(0, M{E t ) - 1), (19) 

where the multiplicity M(Ej) is given by a fit to the experimental data in 
Ref. @ 

M(Ei) = 0.42 + 4.67 lO -4 ^- 1.29 10~ 8 Ef, (20) 



and Ei is given in keV. The motivation of Eq. ( |18| ) is that during the extraction 
method of primary 7 spectra of Ref. |uq] the second and higher generation 7 
ray spectrum, which has of the order M 2 T counts, is subtracted from the total 
unfolded 7 ray spectrum, which has of the order M\ T counts. The errors of 
these spectra are roughly the square root of the number of counts. If we assume 
that these errors are independent from each other, the primary 7 spectra has an 
error of roughly y/ (M\ + M 2 )T. The factor 2 in Eq. ([18]) is due to the unfolding 
procedure and is quite uncertain. We assume this factor to be roughly equal 
the ratio of the solid angle covered by the CACTUS detector array of some 15% 
[ [17| to its photopeak efficiency of some 7% at 1.3 MeV |lq] . We have however 
to apply a couple of minor corrections to Eq. (|1J) . 

Firstly, the first generation method (l5| exhibits some methodical problems 
at low excitation energies. The basic assumption behind this method is that the 
7 decay properties of an excited state is unaffected by its formation mechanism 
e.g. direct population by a nuclear reaction, or population by a nuclear reaction 
followed by one or several 7 rays. This assumption is not completely valid at low 
excitation energies, where thermalization time might compete with the half life 
of the excited state and the reactions used exhibit a more direct than compound 
character. This and some experimental problems like ADC threshold walk and 
bad timing properties of low energetic 7 rays, all described in Ref. |18| , oblige 
us to exclude 7 rays below 1 MeV from the primary 7 spectra. For low energetic 
7 rays above 1 MeV, we increase the error bars by the following rule. For each 
excitation energy bin Ei, we identify the channel with the maximum number 
of counts ch max (this occurs typically between 2 and 3 MeV of 7 energy) . This 



is also the channel with the highest error err max , following Eq. (ttq). We then 
replace the errors of the channels ch below ch max by 

v max _ fa \ 

i + i.O' V ■ (21) 

This formula cannot be motivated by some simple hand waving arguments. We 
feel however, after inspecting several primary 7 matrices, that we estimate the 
systematic error of these spectra quite accurate. 

Secondly, the unfolding procedure [p.4| exhibits some methodical problems 
at high 7 energies. Since the ratio of the photopeak efficiency to the solid 
angle covered by the CACTUS detector array drops for higher 7 energies, the 
counts at these energies are multiplied with significant factors in the unfolding 
procedure. Some channels might nevertheless turn out to contain almost zero 
counts, giving differences in counts between two neighboring channels by two 
orders of magnitude. Since the errors are estimated as proportional to the 
square root of the number of counts, the estimated errors of these channels do 
not reflect their statistical significance. In order to obtain comparable errors 
to neighboring channels we check the errors within one excitation energy bin 
from the 7 energy of ~4 MeV and upwards. If the error drops by more than a 
factor 2, when going from one channel to the next higher one, we set the error 
of the higher channel equal to 50% of the error of the previous one. Also this 
rule cannot be motivated by a simple argumentation. It affects, however usually 
only a very small percentage of channels, and an inspection of several primary 
7 spectra gives us confidence in our error estimation. 

It is now very tedious to perform error propagation calculation through the 
extraction procedure. We therefore decided to apply a simulation technique to 
obtain reliable errors of the g and F functions. For this reason, we add statistical 
fluctuations to the primary 7 matrix. For every channel in the primary 7 matrix, 
we choose a random number r between zero and one. We then calculate x 
according to 

r =7-k/> ( -^>* <22 » 

where a is the number of counts and a the error of this channel. By replacing 
the number of counts a with x, we add a statistical fluctuation to this channel. 
This is done for all channels of the primary 7 matrix, and new g( s > and F^ s ' 
functions are extracted, containing statistical fluctuations. This procedure is 
repeated 100 times, which gives reasonable statistics. The errors in g and F are 
then calculated by 
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2.2.5 Normalizing the level density to other experimental data 

As pointed out above, all solutions of Eq. (g) can be generated from one arbitrary 
solution by the transformation given by Eq. (||) . It is of course discouraging that 
an infinite number of equally good solutions exists, however by comparing to 
known data, we will be able to pick out the most physical one. 

At low excitation energies up to typically 2 MeV for even even nuclei, we can 
compare the extracted level density to the number of known levels per excitation 
energy bin (for a comprehensive compilation of all known levels in nuclei see e.g. 
Ref. p9[). At the neutron binding energy, we can deduce the level density for 
many nuclei from available neutron resonance spacing data. The starting point 
is Eqs. (4) and (5) of Ref. @ 

0F ex P 2VaT7 (2 J + 1) cxp(-(J + l/2) 2 /2a 2 ) 

Q{U ' J) = TT^WTT ^^ ( 25 ) 

itt\ ^ exp2\/a?7 1 . . 

Q{U) = rT^WTTTI^' (6) 

where g(U, J) is the level density for both parities and for a given spin J, and 
g{U) is the level density for all spins and parities; a is the spin dependence 
parameter and a the level density parameter. Assuming that / is the spin of 
the target nucleus in a neutron resonance experiment, the neutron resonance 
spacing D can be written as 

^ = \(g(U n , J = I+ 1/2) + g(U n , J = I- 1/2)), (27) 

since all levels are accessible in neutron resonance experiments, and we assume, 
that both parities contribute equally to the level density at the neutron binding 
energy represented by U„. Combining Eqs. (P5|), ( p6|) and (|27|), one can calculate 
the total level density at the neutron binding energy 

2cr 2 1 

Q\U n ) - -jj- (/ + i) exp(-(7 + l) 2 /2a 2 ) + 7exp(-/ 2 /2 ( T 2 ) ' (28) 

where a 2 is calculated by combining Eqs. (9) and (11) of Ref. [g0| i.e. 

a 2 = 0.0888y/aU^A 2/3 , (29) 

and A is the mass number of the nucleus. It is assumed that a 2 has an error of 
~10% due to shell effects |2Q]. One should also point out, that U n is given by 
U n = B n — P, where B n is the neutron binding energy and P the pairing energy 



which can be found in Tabic III of Rcf. |20 for many nuclei. Unfortunately, 
we cannot compare the calculated level density at the neutron binding energy 
directly with our extracted level density, since due to the omission of 7 rays 
below 1 MeV, the g function can only be extracted up to 1 MeV below the 
neutron binding energy. We will however extrapolate the extracted g function 



with a Fermi gas level density, obtained by combining Eqs. fl26|) and (|2 

1 exp2v / aJ7 

6{) = 12V0l776^/3 a l/2£/3/2 ■ ( 3 °) 

This is done by adjusting the parameters A and a of the transformation given 
by Eq. (||) such, that the data fit the level density formula of Eq. (|o|) in an 
excitation energy interval between 3 and 1 MeV below B n , where in most cases 
all parameters of Eq. (pO) can be taken from Tables II and III of Ref . p3| . This 
semi experimental level density spanning from MeV up to B n is then again 
transformed according to Eq. (||) such, that it fits the number of known levels 
up to ^2 MeV and ~1 MeV for even even and odd even nuclei respectively 
and simultaneously the level density deduced from neutron resonance spacing 
at B n . We have to point out however, that after the fit to known data, the 
extrapolation does not have the functional form of Eq. (|3C|) anymore, due to the 
transformation given by Eq. (|3j) applied to the semi experimental level density. 
Therefore, if necessary, a new extrapolation of the experimental data must be 
performed. 

We have successfully implemented the new extraction method in a Fortran 
77 computer code called RHOSIGCHI jD|. The computer code was compiled 
under a Solaris 2.5.1 operating system running on a Dual UltraSPARC station 
with 200 MHz CPU. The execution time of one extraction is in the order of 
10-20 s. The computer code has ^1200 programming lines, excluding special 
library in and output routines. 



3 Applications to spectra 

3.1 Testing the method on theoretical spectra 

The method has been tested on a theoretically calculated primary 7 matrix. 
The theoretical primary 7 matrix was obtained by simply multiplying a level 
density g to a 7 energy dependent factor F according to Eq. (0). The level 
density was given by a backshifted Fermi gas formula 

g(U) = CU~ 3/2 exp(2v / o[7) (31) 

with U — Ef — P. Below the minimum at U = 9/4a a constant level density 
was used. The 7 energy dependent factor was chosen as 

F(£ 7 ) = C^- 2 . (32) 



In addition, a "fine structure" was imposed on both functions, by scaling several 
~1 MeV broad intervals with factors around 1.5-5. Both model functions are 
shown in the upper half of Fig. O. We extracted the g and F functions from the 
theoretical primary 7 matrix using the excitation energy interval of 4 to 8 MeV 
and excluding all 7 rays below 1 MeV. In the lower panel of Fig. 13 we show 
the ratio of the extracted functions to the theoretical functions. After adjusting 
the extracted quantities with the transformation given by Eq. (g) , we can state 
that the deviation from the input functions is smaller than one per thousand in 
the covered energy range of both functions. Tests of the old extraction method 
showed deviations of the order of 10% to 100% |2]J| . We therefore consider the 
new extraction method to be much more reliable. 

3.2 Testing the method on 172 Yb spectra 

We have tested the method on several experimental primary 7 spectra. We 
will in the following discuss a typical example; the 173 Yb( 3 He,a) 172 Yb reaction. 
The experiment was carried out at the Oslo Cyclotron Laboratory (OCL) at 
the University of Oslo, using a MC35 cyclotron with a beam energy of 45 MeV 
and a beam intensity of typically 1 nA. The experiment was running for two 
weeks. The target was consisting of a self supporting, isotopically enriched (92% 
173 Yb) metal foil of 2.0 mg/cm 2 thickness. Particle identification and energy 
measurements were performed by a ring of 8 Si (Li) particle telescopes at 45° 
with respect to the beam axis. The 7 rays were detected by an array of 28 
5" x 5" Nal(Tl) detectors (CACTUS). More experimental details can be found 
in fllTf . The raw data are unfolded, using measured response functions of the 
CACTUS detector array |1J]. After unfolding, a subtraction method is applied 
to the particle 7 matrix in order to extract the first generation 7 matrix fll5f . 
This primary 7 matrix is taken as the starting point for the extraction method 
presented here. 

In Fig. EL we show the normalized, experimental primary 7 spectra at ten 
different excitation energy bins (data points). The errors of the data points 
are estimated as explained above. The g and F functions were extracted from 
the excitation energy interval 4-8 MeV, excluding all 7 energies smaller than 
1 MeV. The lines are the calculated primary 7 spectra, obtained by multiplying 
the extracted level density g and the 7 energy dependent factor F according to 
Eq. (ph. One can see, that the lines follow the data points very well. It can 
also be seen that the errors of the data points are estimated reasonably giving a 
reduced % 2 of ^0.4. Figure is a beautiful example for the claim, that primary 
7 spectra can be factorized according to the Axel Brink hypothesis |Ll|, [l2j . 

Figure |5j shows how the parameters a and A of the transformation given by 
Eq. (||) can be determined in the case of the 173 Yb( 3 He,a) 172 Yb reaction. The 
extracted g function (data points) is compared to the number of known levels 
[|19| per excitation energy bin (histogram) and to the level density at the neutron 
binding energy, calculated from neutron resonance spacing data |22| (data point 
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in insert). The line in the insert is the extrapolation of the g function up to B n 



according to Section 2.2.5 



In the following, the extracted g and F functions are discussed. Both func- 
tions were already published before, using the old extraction method and some 
fine structure discussed below, could already be seen in the previous publication 
M, Figure g shows the level density g and the relative level density, which is 
the level density, divided by an exponential fit to the data between the arrows. 
The parameters of the fit function 

gm = CeME/T) (33) 

are shown in the lower panel of the figure. In the relative level density one 
can see a small bump emerging at ^2.7 MeV probably due to the quenching 
of pairing correlations |l|, ||. One can also see very nicely the onset of strong 
pairing correlations at 1.0-1.5 MeV of excitation energy. 

In Fig. u\ the 7 energy dependent factor is shown (upper panel). On the 
lower panel, the same data are given, divided by a fit function of the form 

F Rt = CE^. (34) 

This function can be used as a parameterization of 

F(E 1 )=E* x+1 a(E J ), (35) 

where cr(E-y) is the 7 strength function and A is the multipolarity of the 7 tran- 
sition. The fit to the data was performed between the arrows, the fit parameter 
n is given in the lower panel. Since other experimental data is very sparse, we 
did not scale F in order to obtain absolute units. However, the extracted fit 
parameter n is in good agreement with expectations from the tail of a GDR 
strength function at low 7 energies B. In the lower panel a merely significant 
bump at ~3.4 MeV is visible, which we interpret as the Pigmy resonance. 

4 Conclusions 

In this work we have presented for the first time a reliable and convergent 
method to extract consistently and simultaneously the level density g and the 7 
energy dependent function F from primary 7 spectra. The new method, based 
on a least square fit, has been carefully tested on simulated 7 spectra as well as 
on experimental data. In order to normalize the data, we count known discrete 
levels in the vicinity of the ground state and use the level spacing known from 
neutron resonances at the neutron binding energy. 

Compared to the previous projection method H, the least square fit method 
gives the following advantages: The iteration converges mathematically. The 
reproduction of the input level densities and 7 strength functions in simulations 
is much better (almost exact). No tuning of the initial trial function is necessary 
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to obtain a reasonable scaled level density, but the newly derived transformation 
properties of the solution enable the user to normalize the extracted quantities 
with known data. The reduced % 2 is estimated reasonably. The errors of the 
extracted quantities are estimated by statistical simulations. 

We have used the new method to reanalyze previously published data and 
for the analysis of more recent data. Especially the ability to extract absolute 
values of the level density g enables us to perform several new applications 

null- 
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A Proof of Eq. (g) 

The functional form of Eq. (g) opens for a manifold of solutions. If one solution 
of Eq. (g) is found, one can generally construct all possible solutions by the 
following transformation 

§(Ei-E 7 ) = QiEi-EjgiEi-Ej, (36) 

F(E 7 ) = F(E 7 )f(E 7 ). 

The two functions g and / have to fulfill certain conditions, since the set of 
functions g and F are supposed to form a solution of Eq. (0) i.e. 

r ,p FjE^) g{E t - Ry) FjKy) ~g{E t - E 7 ) 

(37) 
Inserting Eq. (B6h one can easily deduce 



Ei 

KEjgtEi-R,) Y, F{E' 1 )g{E l -E' 1 )= (38) 

Ei 

2 f{E' 1 )g{E l ~E' 1 )F{E^)g{E l -E' 1 ). 



Since the right side is independent of E 7 , also the left side must be independent 
of Ely, thus the product of / and g must be a function of Ei only yielding 

f(Ejg(E i -E^) = h(E i ). (39) 
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This condition must of course hold for the case Ei = E 7 . Using the short hand 
notation g(0) = A, one obtains 

Af{E 1 ) = h{E 1 ). (40) 

Inserting this result in Eq. (|38|), one gets 

f(E 1 )g{E t -E 1 ) = Af{E l ). (41) 

Analogously, the condition must hold for the case E 1 — 0, and with /(0) = B, 
one obtains 

Bg(Ei)=Af(Ei). (42) 



Inserting this result in Eq. (|41|), one finally gets 

g(E 1 )g{E l -E 1 )=Ag{E l ). (43) 

We will now show, that the only solution of Eq. (J43|) is an exponential function. 
This proof will involve the limit of Eq. (J43|) for small E~ r However, since g is a 
function of only one variable and the variable Ei is unrestricted in the proof, it 
will be valid for all arguments of g. 

By expanding g in Taylor series up to the first order in E 11 one obtains 

[A + g'(0) EA [giEi) - g< ' {E, b ) EL,] = Agfa). (44) 

Neglecting second order terms in E 7 and dividing by E 7 one gets 

g'(0)g(E i )=Ag'(E i ). (45) 

Defining a — A/g'(0), this differential equation is solved by 

g(E t ) = Ae aE >. (46) 

Using Eq. ((42|), we can easily deduce / to be 

f(E t ) = Be aE \ (47) 

Thus, we have proven the transformation given by Eq. (H) to be the most general 
way to construct all solutions of Eq. (0) from one arbitrary solution. 
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iteration 


P(%) 


number of iterations 


max. variation 


1-5 
6-12 
13-21 
22-30 
31-50 


20 
10 

5 
2.5 

1 


5 
7 
9 
9 
20 


1.2 5 « 2.49 

l.l 7 w 1.95 

1.05 9 w 1.55 

1.025 9 « 1.25 

1.01 20 w 1.22 






E = so 


J]« 11-46 



Table 1: P as a function of the number of iterations. We have in the actual 
computer code implemented 50 iterations where P is decreasing gradually from 
20% to 1%. The maximum variation from the th order estimate any point can 
get is approximately a factor 11. 
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Figure 1: The 7 decay probability from an initial excitation energy Ei in the 
statistical region is proportional to the level density at the final excitation energy 
Ef and the 7 strength function at the 7 energy E 7 = Ei — Ef . 
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Figure 2: Performing the sum of Eqs. ( |lC)| ) and (^lj) along the thick lines. The 
shaded areas are usually excluded when extracting real data due to methodical 
problems (see text). 
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Figure 3: Theoretical level density g and 7 energy dependent factor F used to 
calculate a primary 7 matrix (upper half) . Ratio of extracted to theoretical g 
and F functions (lower half). 
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Figure 4: Normalized experimental primary 7 spectra with estimated errors 
(data points) and calculated primary 7 spectra from the extracted level density 
q and 7 energy dependent factor F according to Eq. (H) (lines). 
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Figure 5: Picking out the most physical solution of Eq. (J2|) by comparing the 
extracted level density g (data points) to the number of known levels per exci- 
tation energy bin (histogram) between the arrows, and the extrapolation of g 
up to B n (line in insert, see text) to the neutron resonance spacing data (data 
point in insert, see text). 
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Figure 6: Extracted normalized level density in 172 Yb (upper panel), and the 
same data divided by an exponential fit (lower panel, see text). 
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Figure 7: Extracted 7 energy dependent factor F in 172 Yb (upper panel), and 
the same data divided by a fit function (lower panel, see text). 
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